Cold dark matter identification: Diurnal modulation reexamined 
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We report on new estimates of the modulation expected in semiconductor detectors due to eclipsing 
of dark matter particles in the Earth. We reevaluate the theoretical modulation significances and 
discuss the differences found with previous calculations. We find that a significantly larger statistics 
than previously estimated is needed to achieve the same level of sensitivity in the modulated signal. 
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I. INTRODUCTION 

For the past ten years or so semiconductor detectors 
(germanium and silicon) have been used to place bounds 
on masses and cross sections of cold dark matter candi- 
dates iB-01 . In these devices massive particles in the halo 
of our Galaxy, with interactions weaker than electromag- 
netic, can be detected by measuring the recoil energy 
produced by their elastic scattering off the detector nu- 
clei. 

To discriminate the signature of these weakly inter- 
acting massive particles (WIMP's) from the background, 
the first idea advanced p] was to look for annual mod- 
ulations of the detection rates originating in the orbital 
motion of the Earth around the Sun. Due to this, the 
average WIMP velocity relative to the Earth does not 
remain constant in time, giving rise to a variation, be- 
tween June and December, of the order of 4% to 6% [g. 

In 1992, Collar and Avignone Q put forward the idea 
of looking for daily — as opposed to annual — modulations 
in counting rates and energy spectra. The argument was 
based on the fact that the Earth, during its daily rota- 
tion, could act as a shield in front of the detector eclipsing 
the dark matter wind traversing the planet. For can- 
didates with scalar or vector interactions (axial-vector 
interactions with the Earth constituents are small) this 
would modify fluxes and velocity distributions, thus af- 
fecting counting rates. The theoretical estimates of this 
effect were obtained by Monte Carlo calculations (MC) 
by these authors [R and indicated that for a suitably 
located detector the modulation could be of the order 
of 2 to 10% for the range of masses for which semicon- 
ductor detectors are more sensitive (10 GeV to 10 TeV) 
and for couplings of the order of the weak coupling con- 
stant H. Either type of search, annual or diurnal, re- 
quires the signal to be collected in the form of "event by 
event;" namely, each event must be time tagged. 

Because of its shorter period, daily modulation has 
the advantage of placing less stringent requirements on 
the stability of the detector and its associated electron- 
ics. The first results ever published of a search for daily 
modulation appeared some time ago ]8|-p0|| and were per- 



formed with a Ge detector in the Canfranc tunnel. The 
experiment, run by the USC/PNNL/UZ (University of 
South Carolina, Pacific Northwest National Laboratories, 
Zaragoza) collaboration was located, however, in a geo- 
graphical location which did not favor this sort of search. 
To improve this, a new site was sought and eventually 
found in Sierra Grande, Argentina, where a laboratory 
was set up and a 1.033 kg Ge detector installed in mid 
1994. The TANDAR group, with support of the rest of 
the collaboration, has been running the experiment since 
then and the first analyses of the data were published 
elsewhere |ll| . 

To get acquainted with the details involved in the cal- 
culation and to carry out an independent check of the 
previous theoretical estimates, a new MC code was de- 
veloped [Q. The new results indicate that, for masses 
larger than ~ 1 TeV (and couplings as large as 10 times 
the weak coupling constant) the modulation of the signal 
is predicted to be less than 1%. This is at variance with 
the original calculations and implies that a larger amount 
of data than previously estimated (more than 20 kg yr for 
a modulation of ~ 0.5% as opposed to 2 kg yr) is needed 
for the sensitivity of the diurnal modulation method to 
exceed that of the conventional signal-to- noise [QJ|] . 

The aim of this paper is to report on new estimates, to 
evaluate the expected modulation significances, and to 
discuss the source of the discrepancy with respect to the 
original calculation. The next section describes the input 
to the calculation and gives a general idea of its contents. 
In the last section we give an account of the results and 
we discuss the differences with previous reports. Finally, 
we present some conclusions. 



II. EXPECTED DAILY MODULATION 

For an incoming WIMP of mass m^ and velocity v 
(/3 = v/c) the total counting rate in the recoil-energy 
interval T, T + dT of the Ge nuclei, is given by 



dN 
dTdt 



(T) = Na 



Phalo 



\'{v)v^{v,T)dv, (2.1) 
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FIG. 1. Relevant directions in the scattering of WIMP's by 
the Earth as it travels around the galactic center with velocity 
Ve- The angle Q defines isodetection rings and a zenith angle. 

where Na is the total number of atoms in the detec- 
tor, phaio/fnx i^ *h^ number density of WIMP's {puaio = 
0.3 GeV cm~^) and g'{v) is the Maxwell-Boltzmann ve- 
locity distribution of the WIMP's for an observer on the 
Earth: namely [O, 



VTT V 2xri 



dv. 
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In the previous equation, 
the Earth with respect to 



Ve is the velocity of 
the galactic rest frame 



{Ve ^ 230 km s~^), ct^, is the velocity dispersion 
of the WIMP's in the galactic halo [a^ is typically 
270 km s^^ pSl), and the distribution (2.2) is normal- 
ized to one. In Eq. (2.1), Vmax is the maximum velocity 



FIG. 3. Difi'erence between velocity distributions of 
WIMP's of mass Mx = 50 GeV for several an- 
gles 6 and the original Maxwell-Boltzmann distribution, 
A = PMc{v,0)dv — PMB{v)dv. Here an interaction strength 
g = IQ was assumed. 



of the WIMP's (w„ 



Ve 



where 



the 



escape velocity from our galaxy, 570 km s"^) and Wmm is 
the minimum velocity of a WIMP necessary to contribute 
to a particular energy of the recoil spectrum, Tmm- If M 
is the reduced mass of the WIMP-nucleus system and 
TUN is the mass of the recoiling nucleus, then 
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In Eq. (2A) the differential cross section for the scat- 
tering of a WIMP off a spin zero (J — 0) nucleus in the 
detector assuming the exchange of Zg-like vector boson 
(G^ = gG\] Gf = 1.16639 x 10"^ GeV-2, the Fermi 
coupling constant) is of the form [Q 
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FIG. 2. Average variation of 9{t) vs LST at Sierra Grande 
throughout the year of 1994. The coordinates of Sierra 
Grande are 41° 40' S, 65° 23' W. 
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FIG. 4. Average velocity of WIMP's of mass Mx=50 GeV 
vs. 6. Three cases g = ?>, g = IQ, and g = 100 are compared to 
the average velocity of WIMP's going through a "transparent" 
Earth. 
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FIG. 5. Energy spectra for several masses at two different 
angles, 9 = 90° (day) and 9 — 170° (night). An interaction 
strength g — 10 was assumed. 
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where R — 1.2yl^/'^ (fm) is the nuclear radius, 9^ is the 
weak mixing angle, and the exponential factor in Eq. 
(2.5) takes into account the loss of coherence of the in- 
teraction. The calculated recoil spectrum must also be 
convoluted with a function accounting for the efficiency 
of the recoiling process generating an ionization signal. 
In our case this function was taken from Ref. |l| . 

Estimations of a diurnal modulation involve consider- 
ing the relation between the rotational axis of the Earth 
and its velocity (see Fig. |^). The velocity of the Earth, 
Ve, through the Galaxy defines an axis of symmetry 
around which, the dispersion of WIMP's has azimuthal 
symmetry. The position vector of a given detector (with 
the origin in the center of the planet) and Ve define an 
angle, 6, that plays the role of a zenith angle. Thus, 
at = 0° the detector is maximally exposed to the 
"WIMP's wind," whereas dX — 180° this exposure is 
diminished by the bulk of the Earth's mass. The angle 
6 also defines isodetection rings, or rings of constant flux 
due to the symmetry of the problem. To each angle 9 we 
can associate a time in some system. Particularly, the 
relation between the local sidereal time (LST) and the 
isodetection angles 9 crossed by the Sierra Grande detec- 
tor during 1994 is shown in Fig. 0. A detectable modula- 
tion effect can be quantitatively predicted provided there 
is a significant change with 9 in the total counting rate 
and/or spectral shape. 

The purpose of the MC code is to track all the parti- 
cles and calculate their number density per unit area and 
their velocity distributions at the surface of the Earth, 
for all isodetection rings. The parameters and model of 
Earth are similar to those used in Ref. 0. A typical 
outcome of the Monte Carlo can be seen in Fig. g| which 
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FIG. 6. Expected rate of events between Emin and 50 keV 
for WIMP's of different masses as a function of local side- 
real time (LST) at Sierra Grande. An interaction strength of 
g = 100 was assumed to make the diurnal modulation effect 
noticeable. 

shows how the velocity distributions, of WIMP's of 
50 GeV mass and 5 = 10 vary for different angles. 
The parameter A depicted in the figure, is defined as 
the difference between the velocity distribution of the 
WIMP's for a given angle 9 and the origin al M axwell- 
Boltzmann distribution as given by Eq. ( |2.2| ), [A — 
Pmc{v, 9)dv — PMB{v)dv]. The distributions are all nor- 
malized to unity, since the number of particles crossing 
the area element dA = 2TTsm9d9 are not the same. Note 
that the elastic scattering with nuclei of the Earth has 
the effect of depleting the region of particles with high ve- 
locities and increasing the number of particles with lower 
speeds. 

A parameter that summarizes the properties of the ve- 
locity distributions is the average velocity distribution 
of incoming (i.e., penetrating the Earth) plus outgoing 
(exiting the Earth) WIMP's, <v>{9). Figure shows 
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FIG. 7. Percentage of modulation, &, as a function of 
WIMP mass for different interaction strengths. 



TABLE I. Predicted diurnal modulation, 5, for several coupling constants. Rmax and Rmin are the total counts expected in 
the energy interval from Emin to 50 keV at maximum and minimum respectively, divided by the energy interval. 
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at max. 
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(keV) 


(counts/keV kg d) 


(counts/keV kg d) 


(%) 


(GeV) 


(h) 


(h) 










ff = 3 


20 


17.5 


9.5 





7.69 


7.51 


2.30 


30 


20.5 


9.5 





10.16 


9.84 


3.12 


50 


22.5 


9.5 





10.55 


10.14 


3.88 


70 


20.5 


9.5 


2.5 


6.95 


6.68 


3.83 


100 


20.5 


8.5 
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5.62 


5.43 


3.43 


200 


20.5 


8.5 


10 


1.73 


1.68 


2.75 


300 


19.5 


8.5 


15 


0.83 


0.82 


1.57 


ff = 10 
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17.5 
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6.04 


5.82 


3.74 


20 
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7.52 


30 
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9.83 


50 
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9.5 
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10.74 


70 


20.5 


9.5 


2.5 


23.10 


20.47 


11.6 


100 


20.5 


9.5 
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18.70 
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10.2 


200 


20.5 


9.5 


6 


8.58 


8.07 


5.9 


300 


20.5 


9.5 


10 


4.28 


4.10 


4.4 


400 


22.5 


9.5 


12 


2.84 


2.74 


3.34 


500 


22.5 


9.5 


15 


1.83 


1.78 


2.83 


1000 


17.5 


11.5 


22 


0.58 


0.57 


1.62 


ff = 100 


10 


19.5 


9.5 





60.94 


44.24 


27.40 


20 


19.5 


8.5 





257.30 


151.04 


41.30 


30 


19.5 


8.5 





335.08 


183.08 


45.18 


40 


19.5 


9.5 





351.34 


188.54 


46.34 


50 


19.5 


8.5 





344.52 


186.06 


45.99 


70 


20.5 


8.5 





308.92 


173.04 


43.99 


100 


20.5 


8.5 





255.74 


152.88 


40.22 


200 


20.5 


9.5 





153.80 


105.88 


31.16 


500 


20.5 


9.5 


5 


43.06 


33.51 


22.20 


700 


20.5 


8.5 


5 


25.15 


20.64 


17.93 


1000 


20.5 


8.5 


10 


14.53 


12.71 


12.56 


2000 


20.5 


10.5 


15 


5.04 


4.74 


5.90 


3000 


19.5 


9.5 


18 


2.72 


2.61 


3.80 



the average velocity of the WIMP's as a function of to- 
gether with the average velocities for the case of a trans- 
parent Earth. Note that for all angles, the average values 
of the distributions are lower than in the transparent case 
meaning that, due to the interaction, the particles lose 
velocity to some degree. In particular, those distributions 
when 9 approaches 180° would yield the largest effect. 

Figures |^ and ^ indicate that the passage of the 
WIMP's through the Earth alters their velocity distri- 
butions. This phenomenon and the net change in the 
number density of particles for rear angles modify the 
differential rate with 9 and are therefore responsible for 
the effect of diurnal modulation. 

Figure illustrates the changes produced in the differ- 
ential energy spectra. In general, low masses are ineffi- 
cient in depositing energy in the detector because their 
momentum transfer to recoiling nuclei is smaller, therefo- 



re their signal decreases rapidly with energy. Signals at 
high energies are only achievable by high-speed WIMP's, 
but WIMPS with such velocities (~ 800 km s"^) are 
scarce. High masses, on the other hand, can deposit 
larger amounts of energy. But, their number density de- 
creases as the mass of the candidates increases because 
of the constant halo density that galactic dynamics re- 
quires. This results in an overall signal decrease with 
WIMP mass at high masses and a different slope on the 
graph. 

The rather small variations in the velocity distribu- 
tions &i 9 = 90° (day) and 6 — 170° (night) are responsi- 
ble for the two sets of predicted energy spectra in Fig. g 
(solid and dashed lines, respectively). Generally speak- 
ing, those WIMP's scattering the Earth at rear angles 
have lower relative speeds and therefore their signal de- 
creases. However, this is only true for a specfic range of 



masses. For masses 100 GeV and larger, the variations in 
the velocity distributions change the slopes of the spec- 
tra generating crossing points at definite energies. These 
crossings energies, Ecross, occur at 1 keV for a mass of 
100 GeV, 11 keV for 400 GeV, 20 keV for 1 TeV, and 
22 keV for a 2 TeV WIMP mass. Determining the min- 
imum energy Emin > Ecross and integrating the energy 
spectra between Emin and an upper hmit (say 50 keV) 
we optimize the predicted signal to obtain the maximum 
variation in the detection rate above Emin ■ 

These new calculations clearly show the energy Ecross 
below which scattering in the Earth increases the count- 
ing rate of WIMP's, and above which the counting rate 
is decreased. This gives a new and more sensitive tech- 
nique to analyze experimental data for diurnal modula- 
tions than just using the portion of the spectrum above 
E 

The variation of the differential rate with 9 is all we 
need to predict the WIMP's signal. By combining this re- 
sult with our knowledge of as a function of time (Fig. |^) 
one is able to look for the characteristic signature in our 
detector due to WIMP's. 



III. RESULTS AND CONCLUSIONS 

The total expected number of counts per keV and kg h 
in the range E^in to 50 keV is shown in Fig. g. In the 
graph, the cross sections were enhanced a factor 100 with 
respect to G^ in order to make the diurnal modulation 
effect visible. From the figure, it can be clearly seen that 
the modulation effect is larger for masses close to 50 GeV, 
whereas it is very small for masses greater than 500 GeV. 
This is due to the asymmetry between the masses of the 
Earth nuclei and the WIMP masses. If the mass of the 
latter is many times larger than that of the nuclei in the 
Earth, the WIMP continues its trajectory almost unaf- 
fected by the interaction and no modulation is expected. 
For masses close to the average mass of nuclei in the 
Earth, the WIMP may suffer a large deflection and the 
scattering angle can be significantly different from zero. 
In this case, we predict a clear modulation effect. 

To quantify the variations in Fig. |6|, we define the mod- 
ulation amplitude, J, as the ratio of the difference be- 
tween the maximum and minimum of the total predicted 
rates (between Emin - 50 keV) with respect to the max- 
imum rate. 



100 



Rn 



R„ 



Rn 



(3.1) 



The behavior of the modulation amplitude, 6, as a 
function of the WIMP mass for different coupling con- 
stants, g, is shown in Fig. M (the actual values are given 
in Table |) . In the graph only data points for which mod- 
ulation is clearly seen are shown. Table I also shows the 
LST at which the maxima and minima are predicted. 
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FIG. 8. Previous predictions of Ref. M (Fig. 9) in solid 
lines and new values obtained after the code was modified 
(dashed lines). 

We notice that for the coupling constants considered 
here, the greater the mass of the candidate, the smaller 
is the modulation expected. This decrease of the modu- 
lation amplitude moves towards the lower masses as we 
decrease the coupling constant. 

The results of the present calculations are in disagree- 
ment with those of Refs. H,0] where no degradation of the 
modulation amplitude is predicted for increasing masses. 
A cross-check of the calculations revealed that in all pre- 
vious results a transformation of the scattering angle 
from the center of mass to the laboratory frame was 
missing. When the correction is included the modula- 
tion amplitudes decrease for masses larger than '^ 50 
GeV bringing the results into agreement. 

Figure g summarizes the differences in predictions at 
1.5 keV before and after the correction. This figure is 
analogous to Fig. 9 in Ref. 0. Notice that, since some 
of the new values of the relative residual are negative. 
Fig. @ shows their absolute values. At some point, the 
spectrum of the high band is larger than the low band 
and this produces a negative value of the residual. Only 
three interaction strengths are considered in the figure 
just to exhibit the modifications to the original and to 
show the behavior of the relative residuals with the mass 
of the WIMP's. Now the new values show a behavior 
with WIMP's mass similar to Fig. 0. 

In view of these new results, we conclude that a larger 
amount of data than previously estimated is needed for 
the sensitivity of the diurnal modulation method to ex- 
ceed that of the conventional signal-to-noisc. The mini- 
mum modulation amplitude that can be detected or ex- 
cluded decreases, as is the case for any search for a pe- 
riodic modulation, as roughly the square root of the vol- 
ume of data. For collected statistics of ~ 500 kg day, and 
typical background levels of 2-3 counts/keV kg day |1^] 
this translates in a sensitivity oi 6 > 2%, which would 
be produced by particles (see Fig. 0) that are already 



excluded using the conventional signal-to-noise method 
(for Dirac neutrinos, g = 1, masses larger than 26 GeV 
and up to 4.7 TeV have already been ruled out by direct- 
detection experiments ||l5[). With these figures, at least 
20 kg yr of data are required for the diurnal modulation 
method to be sensitive to (5 « 0.5% needed to extend the 
WIMP search into the cosmological interesting region. It 
is worth noticing that though the calculations were con- 
ducted assuming vector interactions between the candi- 
dates and the Earth constituents our results also apply 
to dark matter particles with scalar interactions. 

The other alternative to identify a WIMP candidate 
using a semiconductor detector is to look for annual 
modulation of the signal. Now the average velocity of 
the WIMP's impinging on the Earth is what changes 
throughout one period (a year). Since this variation is 
independent of the WIMP's mass, changes in the rate 
should be present for any candidate. At the present time, 
an exhaustive data analysis of the DEMOS collaboration 
is being carried looking for annual modulations of the 
signal in different parts of the spectrum. 
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